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Abstract 

A chaotic model of spontaneous (without external stimulus) neuron firing 
has been analyzed by mapping the irregular spiking time-series into tele- 
graph signals. In this model the fundamental frequency of chaotic Rossler 
attractor provides (with a period doubling) the strong periodic component 
of the generated irregular signal. The exponentially decaying broad-band 
part of the spectrum of the Rossler attractor has been transformed by the 
threshold firing mechanism into a scaling tale. These results are compared 
with irregular spiking time-series obtained in vitro from a spontaneous activ- 
ity of hippocampal (CA3) singular neurons (rat's brain slice culture). The 
comparison shows good agreement between the model and experimentally 
obtained spectra. 
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1. Introduction 

In order to work together the brain neurons have to make an adjustment 
of their rhythms. The main problem for this adjustment is the very noisy 
environment of the brain neurons. For pure periodic inner clocks this adjust- 
ment would be impossible due to the noise. Nature, however, has another 
option. This option is a chaotic clock. In chaotic attractors certain char- 
acteristic frequencies can be embedded by broad-band spectra, that makes 
them much more stable to the noise perturbations [lj]. We will present em- 
pirical evidence supporting the proposition that the Rossler chaotic attractor 
in combination with an appropriate threshold passage mechanism could be 
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Figure 1: Mapping of a spike train (figure lb) into a telegraph signal (figure la). 



used as a simple model of the spontaneous (without external stimulus) firing 
activity frequently observed in in vitro hippocampal neurons 

The spontaneous firing activity should be more simple and self-consistent 
than a reaction of a neuron to the external stimuli. This can allow an anal- 
ysis of the inner neuron clock in its "free-run" mode. Consideration of the 
most frequently firing neurons can also help in this direction. Besides its 
pure academic significance for studying the neurons' firing the sponta- 
neous activity provides a considerable contribution to network development 
[Bt] , , information processing Q- 10 1 , and behavioral responses U|, 12]. It 
is known that periodic spontaneous bursts of the activity can convey infor- 



mation about sensory stimuli [13j],|l4|]. Spontaneous activity in brain slice 



preparations purely reflects the intrinsic properties of local circuits and indi- 
vidual neurons and hence allows for the investigation of the internal dynamics 



of neuronal networks 15-18 



All types of information, which is received by sensory system, are en- 
coded by nerve cells into sequences of pulses of similar shape (spikes) before 
they are transmitted to the brain. Brain neurons use such sequences as main 
instrument for intercells connection (both for tuning of their chaotic clocks 
and information- wise). The information is reflected in the time intervals be- 
tween successive firings (interspike intervals of the action potential train, see 
Fig. lb [l9j]). There need be no loss of information in principle when con- 



verting from dynamical amplitude information to spike trains [20[ and the 
irregular spike sequences are the foundation of neural information process- 
ing. Although understanding of the origin of interspike intervals irregularity 
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Figure 2: X-component fluctuations of a chaotic solution of the Rossler system Eq. (f ) 
(a = 0.15, 6 = 0.20, c= 10.0). 



has important implications for elucidating the temporal components of the 
neuronal code and for treatment of such mental disorders as depression and 
schizophrenia, the problem is still very far from its solution (see, for instance, 



Ref. 21( and references therein). Motivation to study the hippocampus in 



relation to depression is based on the recently discovered evidences of its 
deep involvement in this mental disorder. The hippocampus is a significant 
part of a brain system responsible for behavioral inhibition and attention, 
spatial memory, and navigation. It is also well known that spatial memory 
and navigation of the rats is closely related to the rhythms of their moving 
activity. On the other hand, the hippocampus of a human who has suffered 
long-term clinical depression can be as much as 20% smaller than the hip- 
pocampus of someone who has never been depressed (22| . One can speculate 
that in the case of depression the chaotic neuron clocks can be broken in a 
significant part of the brain neurons. That can result in certain decoherence 
in different parts of the brain. 

Treatment of the spiky signals by analytical methods also presents a dif- 
ficult problem. In order to apply the Fourier transform method one can use 
different mappings: mapping of the finite spikes into Dirac delta functions 
or into telegraph signals, for instance (cf. Ref. [23]). For a single neuron 
firing the amplitude of the spikes are almost identical to each other and the 
neural information is coded in the length of the interspike intervals and their 
positions on the time axis 



19 , 20 



therefore it is the most direct way to map 
the spike train into a telegraph time signal, which has values -1 from one side 
of a spike and values +1 from another side of the spike with a chosen time- 
scale resolution. An example of such mapping is given in figure 1. While the 
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Figure 3: Z-component fluctuations of a chaotic solution of the Rossler system (the pa- 
rameters are the same as in Fig. 2). 



information coding is here the same as for the corresponding spike train, the 
Fourier transform methods are quite applicable to analysis of the telegraph 



time-series 



24 



2. A spontaneous firing model: the Rossler system 

Nerve cells are surrounded by a membrane that allows some ions to pass 
through while it blocks the passage of other ions. When a neuron is not send- 
ing a signal it is said to be "at rest". At rest neurons exhibit very small con- 
ductance of sodium ions and slightly larger potassium conductance against a 
high concentration of intracellular potassium ions. The resting value of mem- 
brane electrochemical potential - P (the voltage difference across the neural 
membrane) of a neuron is about -70mV. If some event (a stimulus) causes the 
resting potential to move toward OmV and the depolarization reaches about 
-55mV (a "normal" threshold) a neuron will fire an action potential. The 
action potential is an explosive release of charge across plasma membrane 
and its surrounding that is created by a depolarizing current. If the neuron 
does not reach this critical threshold level, then no action potential will fire. 
Also, when the threshold level is reached, an action potential of a fixed size 
will always fire (for any given neuron the size of the action potential is always 
the same). 

Recent reconstructions of a driver of the membrane potential using the 
neuron spike trains indicate the Rossler oscillator as the most probable (and 



simple) candidate (see, for instance, Refs.[25||-[3l|). Figure 2 shows as exam- 
ple the x-component fluctuations of a chaotic solution of the Rossler system 

El 
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Figure 4: Spectrum of the telegraph signal corresponding to the spike train generated by 
the x-component fluctuations overcoming the threshold x — 7. The dashed straight line 
indicates a power law Eq. (3) in the log-log scales. 



dx dy dz 

- = -( y + z>, Tt =^ + ay; Tt =b + xz-cz (1) 

where a, b and c are parameters. At certain values of the parameters a,b 
and c the z-component of the Rossler system is a spiky time series Fig. 3 



(see also Refs. [331] .[34 

It can be shown that the Rossler system and the well known Hindmarsh- 
Rose model [35[ of neurons are subsystems of the same differential model with 
a spiky component [34]. Previously the 'spiky' component of such models was 
interpreted and studied as a simulation of a neuronal output. For the spon- 
taneous neuron firing (without external stimulus), however, we suggest to 
reverse the approach and consider the spiky variable as the main component 
of the electrical input (which naturally should have a 'spiky' character, see 
above) to the neuron under consideration. For any given neuron the height 
of the spikes, which the neuron generates, is about the same. However, the 
heights of the spikes generated by different neurons are different. Also the 
signals coming from different neurons to the neuron under consideration have 
to go through the electrochemical passes with different properties. Therefore, 
the spiky z-time-series (Fig. 3) can naturally represent a multineuron signal, 
which can be considered as a spontaneous input for the neuron under con- 
sideration. If we use the usual interpretation of the x-component as a driver 
of the membrane potential P(x) and the y-component as that taking into 
account the transport of ions across the membrane through the ion chan- 



nels 35|, then addition of the spiky z-component (representing in present 
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Figure 5: Spectrum of the x-component fluctuations shown in Fig. 2. We used the semi- log 
axes in order to indicate exponential decay of the spectrum. 

model the multineuron spontaneous input) in the dynamical equation for x- 
component is similar to the addition of an external input-component to the 
dynamical x-equation in the Hindmarsh-Rose model. Then, the quadratic 
nonlinearity in the third equation of the system Eq. (1) can be interpreted 
as a simple (in the Taylor expansion terms) feedback of the neuron to the 
main component of the neuronal input. This model with the strong nonlinear 
feedback can be relevant to the most active neurons of a spontaneously active 
brain (see below results of an in vitro experiment with a spontaneous brain 
activity). The details of the function P(x) is not significant for the threshold 
firing process, what really matters is that the membrane potential function 
P(x) reaches its firing value when (and only when) its argument x crosses 
certain threshold from below. In this simple model the driving variable x 
may overcome its threshold value (Fig. 2) due to the deterministic (chaotic) 
spontaneous stimulus. Let us consider an output spike signal resulting from 
overcoming a threshold value x — 7, for instance. Fig. 4 shows spectrum of 
the telegraph signal corresponding to the spike train. In order to understand 
what is going on here we show in figure 5 spectrum of the x-component itself. 
The semi-log scales are used in these figures in order to indicate exponen- 
tial decay in the spectra (in the semi-log scales this decay corresponds to a 
straight line): 

E(f) ~ e- ///e (2) 
While the low-frequency peak in the spectrum corresponds to the fundamen- 
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Figure 6: Autocorrelation function for the telegraph signal corresponding to the cell-21 
(800 spikes). The insert shows corresponding spectrum. 



tal frequency, / , of the Rossler chaotic attractor, the rate of the exponen- 
tional decay (the slope of the straight line in Fig. 5) provides us with and 
additional characteristic frequency f e Eq. (2). It should be noted that for a 
wide class of deterministic systems a broad-band spectrum with exponential 



decay is a generic feature of their chaotic solutions Refs. 36]- [39 . 

Thus Rossler chaotic attractor has two clocks: fundamental with fre- 
quency /o and decaying with frequency f e . If one compares Fig. 4 and Fig. 
5 one can see that the fundamental clock survives the threshold crossing 
(with a period doubling, see also Fig. 9). The decaying clock, however, does 
not survive the threshold crossing: the exponential decay in Fig. 5 has been 
transformed into a scaling (power law) decay in Fig. 4 

E(f) ~ r a (3) 

(with a ~ 1.7), which has no characteristic frequency (scale invariance). 



In the above presented model the fundamental frequency of chaotic Rossler 
attractor provides (with a period doubling) the strong periodic component 
of the generated irregular signal. This periodic component can be utilized by 
Nature as a chaotic clock of the spontaneous neuron firing (see next section). 
The exponentially decaying broad-band part of the spectrum of the Rossler 
attractor has been transformed by the threshold firing mechanism into a scal- 
ing tale. To understand mechanism of this transformation is rather difficult 
problem since the mechanism of generation of the exponential spectrum by 
the chaotic systems is itself an unsolved problem (see, for instance, Refs. 
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Figure 7: Spectrum of the telegraph signal corresponding to the cell-21 (800 spikes) in 
log-log scales. The dashed straight line indicates a power law Eq. 4: E(f) ~ f~ 17 - 



36j-[39]). The scaling exponent value '-1.7' is not sensitive to a reasonable 



variation of the threshold value (~ 20%) and even to Gaussian fluctuations 
of the threshold value. Therefore, it is not just a coincidence that the scaling 
law in the Rossler case agrees with results of the in vitro experiment reported 



in the next section (cf. also [24|.|40j]). The power-law spectrum Eq. (3) with 
the exponent a ~ 1.7 is an indication of a strongly non-ergodic system, with 
perennial aging, or maybe that an infinitely aged correlation function is pre- 
dominant due to noise or truncations. Moreover, the in vitro data, which 
we used for comparison (see next Section) seem to fit a renewal hypothesis 



compatible with this scaling exponent, and the theory of a recent Ref. [41 
(see also Ref. (iif). 



3. In vitro spontaneous brain activity 

In order to compare this simple model consideration with the experi- 
mental data we have analyzed spike trains obtained in vitro from a spon- 
taneous activity in CA3 hippocampal slice culture of a Wistar/ST rat (the 
raw data and the detail description of the experiment can be found online 



at http://hippocampus.jp/data and in Refs. [2j,[3|]), In the in vitro exper- 
iment a functional imaging technique with multicell loading of the calcium 
fluorophore was used in order to obtain the spike trains of spontaneously 
active singular neurons (hippocampal pyramidal cells) in the absence of ex- 
ternal input. In this experiment different levels of activity were observed for 



different neurons [2], [3]. We take for our analysis the two most active neu- 



rons (http:/ /hippocampus.jp/data - Data-006, cell-21, with 800 spikes in the 
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Figure 8: As in Fig. 7 but for cell-25 (D-006, 692 spikes). 

time-series; and a cell-25, with 692 spikes). The spike trains were mapped to 
telegraph signals as it is described above. Figure 6 shows the autocorrelation 
function for the telegraph signal corresponding to the cell-21 (800 spikes). 
Insert in the Fig. 6 shows the corresponding spectrum. Both the correlation 
function and the spectrum provide a clear indication of a strong periodic 
component in the irregular signal (the oscillations in the correlation function 
and the peak in the spectrum). The periodic component can be seen at fre- 
quency /o ~ 0.3Hz. Figure 7 shows the spectrum in log-log scales. One can 
see that at high frequencies the spectrum exhibits a scaling behavior Eq. (3) 
(power law: E(f) ~ f~ 1,7 , as indicated by the dashed straight line). The real 
power law can be more pronounced but under the experimental conditions 
individual spikes emitted at firing rates higher than 5Hz were experimentally 
inseparable |2[,[3[. Figure 8 shows spectrum of the telegraph signals corre- 
sponding to the spike train obtained for the cell-25 (D-006, 692 spikes). The 
spectrum is rather similar to the spectrum shown in Fig. 7 ( for cell-21). 
The more broad peak in Fig. 8 can be related to the poorer statistics for the 
cell-25 (number of spikes 692) in comparison with cell-21 (number of spikes 
800). One can compare Figs. 7 and 8 with Fig. 4 to see very good reproduc- 
tion of the main spectral properties. Figure 9 shows also a superposition of 
the autocorrelation functions for the telegraph signals corresponding to the 
cell-21 (the solid line) and to the spike train generated by the the Rossler 
attractor fluctuations overcoming the threshold x = 7 (circles). In order to 
make the autocorrelation functions comparable a rescaling has been made 
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Figure 9: Autocorrelation functions for the telegraph signals corresponding to the cell- 
21 and to the spike train generated by the Rossler attractor fluctuations overcoming the 
threshold x — 7 (circles). In order to make the autocorrelation functions comparable a 
rescaling has been made for the Rossler autocorrelation function. 

for the Rossler autocorrelation function: the scaling coefficient is equal to 
0.024 (the time step size for the Rossler solution was 0.01). As in the case of 
the spectra the comparison of the autocorrelation functions is not sensitive 
to a reasonable variation of the threshold value ( 20%) and even to Gaussian 
fluctuations of the threshold value. 



4. Conclusion 

Simple Rossler system has properties, which allow to use it as an adequate 
model of spontaneous neuron activity (that includes a multineuron input rep- 
resented by the spiky z-component of the Rossler attractor, Fig. 3, and a 
strong nonlinear feedback) . A mapping of the spike time series into telegraph 
signals (that preserves completely intact the 'frequency modulated' informa- 
tion of the spike series) allows to compare spectral properties of the model 
threshold firing (Fig. 4) with analogous firing of the in vitro brain neurons 
spontaneous activity (Figs. 7 and 8). This investigation can be considered as 
an additional confirmation and specification of the already existing evidences 



(Refs. [25]- [31]) that the Rossler system can be a driver of spontaneous neu- 
ron firing. On the other hand, it is known that the spontaneous firing of 
brain neurons accounts for about 80% of the metabolic energy consumed by 



the brain [43]. Therefore the spontaneous neural activity should have a very 
significant neurobiological function (cf. Introduction). 

Spontaneous neural activity with multineuron nonlinear interactions has 
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been studied mainly in neuronal networks (see, for instance, Refs. |2]-[18j). 
In present simple model a multineuron input and a strong nonlinear feedback 
are simulated in the frames of a single self-consistent low- dimensional system 
with a chaotic solution. On the other hand, we do not know any study where 
a fully chaotic system could make anomalous behaviors emerge, so it seems 
that studying threshold passage is enough to have that effect. 

Of course, in vivo neuron signals can be much more complex. The rela- 
tively simple firing of the in vitro spontaneously active hippocampal neurons 
can be used in order to reveal the underlying neuron dynamics. 
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